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Polyelectrolyte-surfactant complex: 
phases of self-assembled structures 
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We study the structure of complexes formed between ionic surfactants (SF) and a single oppositely 
charged polyelectrolyte (PE) chain. For our computer simulation we use the "primitive" electrolyte 
model: while the polyelectrolyte is modeled by a tethered chain of charged hard sphere beads, 
the surfactant molecules consist of a single charged head bead tethered to a tail of tethered hard 
spheres. A hydrophobic attraction between the tail beads is introduced by assuming a Lennard- Jones 
potential outside the hard-sphere diameter. As a function of the strengths of both the electrostatic 
and the hydrophobic interactions, we find the following scenario: switching on and increasing the 
electrostatic forces first leads to a stretching of the PE and then by condensation of SF to the 
formation of a complex. For vanishing hydrophobic forces this complex has the architecture of a 
molecular bottle-brush cylindrically centered around the stretched PE molecule. Upon increasing 
the hydrophobic attraction between the SF tails, a transition occurs inverting this structure to a 
spherical micelle with a neutral core of SF tails and a charged corona of SF heads with the PE 
molecule wrapped around. At intermediate hydrophobicity there is a competition between the 
two structures indicated by a non-monotonic dependence of the shape as function of the Coulomb 
strength, favoring the cylindrical shape for weak and the spherical micellar complex for strong 
interaction. 

PACS numbers: 82.70.-y,61.20.Ja,82.35.Rs,83.80.Qr 



I. INTRODUCTION 



Polyelectrolyte-surfactant mixtures have proven to 
provide the basis for new materials with extraordinary 
properties that make them interesting for a wide range 
of applications 0,0- In particular the Coulomb at- 
traction between polyelectrolyte chains and oppositely 
charged ionic surfactant molecules in solution leads to 
aggregation and complex formation. The resulting com- 
plex differs in its conformational, structural and dynam- 
ical features from that of solution containing only one 
of its constituents, the pure polyelectrolyte or the pure 
surfactant (for reviews see |lj, |3j). Much experimental 
work has been done on these systems and a rich variety 
of different complex structures was revealed as a function 
of chemical nature of the polyelectrol yte , surfactant and 

solvent molecules HHHBHUHElEEEim- 

A theory that could give a systematic microscopic char- 
acterization and prediction of the different complexes on 
the other hand does not exist so far. 

Already for complexes formed by a single polyelec- 
trolyte chain with ionic and hydrophobic surfactants at 
low concentration a theoretic treatment that takes the 
long-ranged Coulomb interaction and the hydrophobic 
interactions fully into account is still to be developed. 
The case without hydrophobicity has recently been sim- 
ulated by the present authors [l6j . For these single poly- 
electrolyte complexes one expects structures either of 
bottle-brush shape for low hydrophobicity or, for strong 



hydrophobicity, possibly spherical micelles that aggre- 
gate together with the polyelectrolyte chain 0, [^j . The- 
ories so far in general make assumptions about the sym- 
metry or the structure of the complex implying a com- 
plete adsorption of surfactant onto the polyelectrolytes 
for rigid [Islllll an d flexible chains, see e.g. [2|| and can- 
not predict the complex structure or its symmetry 
as far as it enters the theory as an input. For the spe- 
cial case of a system of rigid polyelectrolytes that form 
complexes via thermodynamic counterion condensation 
including the effect of added surfactant a theory has re- 
cently been developed by Kuhn, Levin and coworkers 
pll E5 , E3 | . The collapse and partial collapse of semi- 
flexible and flexible polymers in solutions with surfac- 
tants has been treated by a meanfield association the- 
ory by Diamant and Andelman pil |25|. However, this 
theory does not include explicitly the long range of the 
Coulomb interactions. A self consistent field theory with 
a spherically symmetric setup [T^ on the other hand pre- 
dicts the formation of spherically symmetric micelles with 
a neutral hydrophobic core and a surface layer where 
the charged surfactant heads and the oppositely charged 
polyelectrolyte are located. Based on this observation 
it has been proposed to simplify the model by replacing 
the ionic surfactant by large spherical counterions not 
treating the surfactant tails explicitly [2(| . For a system 
of a polyelectrolyte in the presence of small counterions 
computer simulations [53, 0, and theories [3(1 El 
find a swelling and stretching of the polyelectrolyte and 
a subsequent collapse when increasing the Coulomb in- 
teraction. In this scenario a non-monotonic dependence 
of the of the polyelectrolyte 
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FIG. 1: Snapshots of the complex configurations showing 
the dependence on the two interactions: Coulombic (As) and 
hydrophobic (elj) 



conformation as a function of the Coulomb strength oc- 
curs: for weak Coulomb strength the PE stretches while 
in the strong Coulomb regime counterion condensation 
and a collapse to a coil conformation occurs. For large 
counterions representing micelles the polyelectrolyte may 
be wrapped around the charged sphere exhibiting differ- 
ent conformations dep en ding; on the chain flexibilit y |2 fit 

m s h m m ej m m, in m m sHmiiir 

The bottle brush structure that is expected for the 
condensation of charged ionic surfactants onto a poly- 
electrolyte chain may be compared with corresponding 
neutral molecular bottle brushes in case the polarization 
effects in the complex are neglected. Such neutral bottle 
brush or comb polymers are built either by end linking 
the side chains to the backbone chain or by attaching 
them via strong hydrogen bonds 0, HH ■ Such configu- 



rations may be expected to become more rigid for higher 
densities of the side chains. A quantitative description, 
however, is still under debate since different simple scal- 
ing arguments [4^ . IHol l5l| appear to predict different 
quantitative behavior for the sizes of the main and side 
chains depending on the model used and the limit that 
is considered [52| . 

In the present paper we explore by computer simu- 
lation the phases of complex formation of a single poly- 
electrolyte chain with ionic surfactants in the phase space 
parameterized by the strengths of the Coulombic interac- 
tions and the hydrophobic attraction between the surfac- 
tant tails. As in our previous study 0] our simulation 
relies on the so-called "primitive" model [5^ that treats 
the microscopic charges of the polyelectrolyte and the 
surfactant heads explicitly. This is in contrast to Simula- 
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tions on polymer-surfactant complexation which ne glec t 
the long-range nature of the Coulomb interactions |54| . 
The "primitive" model has recently also been applied to 
simulate the complexation of polyelectrolyte chains with 
charged spheres representing spherical micelles [5^ and 
has proven to be reliable in different contexts of polyelec- 
trolyte conformations [H |H El El HHH H| . 

To treat surfactants with neutral or hydrophobic tails 
we model the surfactant by a tethered hard sphere chain 
with a charged head and a Lennard- Jones like effective at- 
traction between the hydrophobic tail monomers allowing 
for the formation of micelles. As a result, we observe dif- 
ferent phases of complexation depending on the strengths 
of the two interactions involved: the Coulomb and the 
hydrophobic interaction. An overview of these phases 
is given in Fig. \T\in terms of simulation snapshots. For 
vanishing hydrophobicity (leftmost column of Fig. ^) the 
mean square end-end distance R 2 of the polyelectrolyte 
is non-monotonic as function of the Coulomb interaction. 
However, in contrast to the standard situation where a 
collapse of the pol yele ctrolyte due to counterion conden- 
sation occurs [27], IH |2j| , such a collapse is blocked by 
the steric interaction of the surfactant tails. As we have 
shown in Ref. 16], for even longer surfactant tails the 
collapse is more sufficiently blocked, while the addition 
of salt weakens the complex as it replaces the surfactant 
in the complex. Increasing the hydrophobicity for van- 
ishing Coulomb interaction on the other hand (bottom 
row of Fig^l leads to the aggregation of micellar clus- 
ters of surfactant molecules that are independent from 
the polyelectrolyte chain. Introducing a finite Coulomb 
strength in this latter situation favors the condensation 
of the surfactant at the polyelectrolyte. Thirdly we may 
start with a bottle brush conformation of the complex 
at high Coulomb but vanishing hydrophobic interactions 
(top left of Fig^). Gradually increasing the hydropho- 
bicity (along the top row of FigQ the tails start to ag- 
gregate and and a transition occurs from the cylindrical 
structure with the tails on the outside to an inverted 
spherical structure with a core that is constituted by the 
surfactant tails. In our present resolution this transition 
appears to be sharp. 

The paper is organized as follows: in section II, we 
describe our model in detail and we present details of 
the simulation procedure in section III. Results are given 
in section IV and we conclude in section V. 



II. THE MODEL 

In our simulation we model all molecules as chains of 
tethered hard spheres (beads) that represent Kuhnian 
segments of the molecule. Microions are modeled by hard 
sphere beads. Charges on these molecules and ions are 
represented by placing a point charges at the center of 
any charged bead. In this primitive model the solvent 
is not treated explicitly and only enters via its dielec- 



tric constant e r neglecting the discreteness of the solvent 
molecules. The microscopic interactions between the par- 
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FIG. 2: Model of the polyelectrolyte - surfactant system 
with salt ions. Charges are indicated by "+" and "-" signs. 

tides are given by the hard sphere excluded volume, the 
long-ranged Coulomb forces, and a short range hydropho- 
bic interaction. We reduce the model parameter space 
by fixing the hard sphere diameter of all beads in the 
system to the same value a. Furthermore, all charged 
beads have the same charge ±qe with e denoting the el- 
ementary charge. While this "minimal" model explicitly 
handles the counter ions and may be used to describe 
strongly charged, flexible polyelectrolytes (2?|], analytical 
calculations usually neg lect the individual character of 
the counter ions [6al6lTl6ll6^|. 

As mentioned above, all molecules are modeled by 
chains of freely jointed, tethered hard sphere beads. Sub- 
sequent beads in the same chain are tethered. The exten- 
sion of all bonds is limited by the tether to a fixed value b. 
For vanishing Coulomb and hydrophobic interactions the 
model is purely entropic (i.e. temperature- independent). 
Our specific parameters are such that the polyelectrolyte 
is represented by a chain of N charged beads each with 
charge — qe. The surfactant molecules on the other hand 
are chains of n tethered beads with one charged head 
bead of charge +qe and n — 1 neutral tail beads that 
may interact due to hydrophobicity. The full model may 
also contain salt which has been treated in Ref.[16(. The 
salt ions are then represented by single charged beads. 
The constituents of the model of polyelectrolytes, surfac- 
tant molecules and microions is schematically shown in 
Figure |21 

In our model, the effective interaction due to hydropho- 
bicity of the surfactant tails is given by a Lennard- Jones 
like potential between the tail monomers. This approach 
has also been used in the context of copolymer micel- 
lization j^, |(55|]. The explicit interaction between any 
two beads i,j of the model at a distance r between their 
centers is given by a potential Vij(r). To denote chain 
connectivity we use an adjacency matrix Ty. For beads 
i,j that are tethered = 1 otherwise TJy = 0. With 
qi G {±1,0} denoting the charge and Ui S {1,0} the 
hydrophobicity number of bead i, the potential Vij(r) is 
given by 
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^i(r) = <M;A B ~ + u lUj e LJ (^) " - (^) 6 ) + F H (a - r) + r«V H (r - 6). (1) 



Here, the dimensionless Coulomb strength Ab is the Bjer- 
rum length measured in bead diameters a. The Bjerrum 
length ctAb = q 2 e 2 /e r k-QT sets the length scale where the 
Coulomb pair interaction is comparable to the thermal 
energy fc^T. The dimensionless hydrophobicity parame- 
ter el j controls the depth of the Lennard- Jones poten- 
tial. Furthermore, the hard core and tether interactions 
are represented by the one dimensional hard wall poten- 
tial 

*<•>-{« ir° ■ p> 

We note that in our notation the product TijVn(r — b) is 
defined to be zero if 2y is zero. 

The key parameters of our model are: the number iV 
of charged monomers of the polyelectrolyte chain, the 
number n — 1 of neutral beads of the surfactant tails, 
the parameters Ab and Slj which measure the strengths 
of the Coulomb and hydrophobic interactions and finally 
the relative salt concentration c s . In the following we 
shall explore the parameter space in particular in the 
two variables Ab and Elj while keeping fixed the par- 
ticle hard-core diameters, the microion charges qi, the 
monomer number TV = 32 and the maximal tether length 
b = f.5cr and in general also the salt concentration c s = 
and the surfactant length n = 5. In our previous work 
[l6j we also investigated the influence of surfactant tail 
variation as well as of finite salt concentration. A brief 
discussion of the corresponding results is also included in 
the following. 



III. SIMULATION TECHNIQUE 

We simulate by standard Monte Carlo (MC) methods 
a system of a single polyelectrolyte chain of N charged 
beads at finite concentration with polyelectrolyte charge 
density p = N/L 3 = O.OOlcr -3 which is kept fixed 
throughout all simulations. The finite concentration is 
taken into account by periodic boundary conditions of 
our finite cubic simulation box of length L and replicated 
images of the charges. Summations over the interactions 
with the periodic images are calculated by the Lekner 
sum method • The Metropolis rates of the MC moves 
are determined by the interactions given by Eqn. 

The relaxation of the system, especially in the fully as- 
sembled bottle brush configurations is very slow. This is 
also known from simulations of conventional bottle brush 
molecules 52]. To ensure sufficient relaxation we per- 
formed ~ 5 ■ 10 6 attempted MC moves per particle at 
each state point. The acceptance ratio is roughly 0.8 in 
situations with an open structure. In the micellar cases 
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however, caging occurs inside the densely packed core 
reducing the acceptance rate down to 0.3. 

In all simulations the global charge of our system van- 
ishes. We performed simulations for a system of a single 
polyelectrolyte chain together with TV = 32 oppositely 
charged surfactant molecules. The case n = 1 of no tail 
beads at all serves as a natural reference for pure counter 
ions. Such a situation was already studied in Ref. |27j 
and implemented in the frames of the present setup in 
Ref. for different Bjerrum lengths. Note, however, 
that chain connectivity was modeled differently in Ref. 
|27| via a finite extension potential. 

IV. RESULTS 

A. Zero hydrophobicity of the surfactant 

In the case of vanishing hydrophobic interactions 
£lj = the surfactant tails interact only due their steric 
repulsion. This special case was treated by the present 
authors in |lq in detail, so we only recall some of the 
basic results. We first discuss the case of no added salt. 
As a reference situation serves the a surfactant without 
tail, i.e. n = 1. The averaged square end-to-end distance 
R 2 of the polyelectrolyte chain is defined via 

R 2 =< (i?i - R N f > (3) 

where < ... > denotes a statistical average and R\ and 
Rn are the actual positions of the two end- monomers of 
the polyelectrolyte. 

R 2 is plotted as a function of Coulomb strength Ab 
in Fig. |3J Here, the collapse scenario found in other 
simulations [27l 123 . 129) and analytical treatments [(57| is 
confirmed. This scenario for simple countcrions is the fol- 
lowing: At zero Bjerrum length As = the interaction 
is purely excluded volume and the chain attains a poly- 
mer coil configuration. Increasing Ab the chain stretches 
due to the repulsion of its charged beads. For low As 
the counterion screening is weak and stretching contin- 
ues until <j\b is of the order of the counterion diameter, 
i.e. Xb — f ■ For larger Xb > 1 the condensation of coun- 
terions induces a shrinking of the PE extension which for 
high Ab may even be become smaller than the neutral 
coil; an effect not yet seen in the range for Ab that we 
simulated. 

Substituting in the former system the counterions by 
charged surfactants with neutral tails the behavior as 
function of the Coulomb strength Ab changes in the fol- 
lowing way: The PE chain stretches for increasing Ab 
with the maximal extension attained for higher values of 
Ab corresponding the the larger effective diameter of the 
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FIG. 3: The mean square end-to-end distance R as a func- 
tion of the Bjerrum length Xb for n = 1(»), 5(D) and 10(a) 
(from below) monomer surfactants for vanishing hydropho- 
bicity parameter el.j — 0. The statistical uncertainty is indi- 
cated by an error bar. 



FIG. 5: The number of clusters N c i in the system as func- 
tion of the Coulomb strength Xb for different values of the 
hydrophobicity parameter el.j = 0(D), 1(B), 2(o), 4(»), 6(A 
),8(a) from above. 
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FIG. 4: The mean square end-to-end distance R 2 for Bjerrum 
length \b — 8a as a function of the relative salt concentration 
c 3 for (from below) n = 1(»), 5(D) and 10(A) monomer sur- 
factants. The statistical uncertainty is indicated by an error 
bar. 



vor of salt ions. This is driven by entropy of mixing and 
corresponds to a chemical equilibrium. The replacement 
is therefore increasing with increasing salt concentration. 
This destroys the structure of the complex as seen from 
the reduction in the end-to-end distance of the polyelec- 
trolyte. 

In all our simulations we observe that the contour 
length L of the polyelectrolyte given by the sum of the 
bond lengths is very stable under all changes of the 
parameters. Varying the Bjerrum length in the range 
Xb = ... 8 we find a corresponding increase from L = 40 
to L = 42, a change that is within the statistical uncer- 
tainty. For this reason also the persistence length L p does 
not give additional information. For a worm-like chain 
namely it is related to the end-end distance R and the 
contour length L by [63 

R 2 = 2LL p - 2L 2 p (l - exp(-i/L p )) (4) 

a relation also confirmed by simulations [27J • 



charged surfactant acting as a counterfoil. When coun- 
terion condensation occurs at higher Xb the surfactant 
and the PE chain form a complex in which the surfac- 
tants aggregate with their heads near the PE chain and 
their tails pointing away from the PE chain in a man- 
ner resembling the structure of a molecular bottle brush. 
The subsequent collapse of the PE chain that is observed 
in the previous scenario now is effectively blocked due to 
the excluded volume interaction of the surfactant tails. 
The internal structure of this complex is discussed be- 
low. For the special case of vanishing hydrophobicity the 
influence of added salt on the end-to-end distance R was 
investigated in Ref.0. Adding N s additional 1 : 1 salt 
pairs to our system we define the relative concentration 
as c s — N s /N where N is the number of charges on the 
polyelectrolyte. The results are shown in Fig. 0] The 
surfactant molecules in the complex are replaced in fa- 



B. Clusters and micelles 

The results described in this and the following sections 
section all refer to a system with a single polyelectrolyte 
chain of 32 charged beads and 32 oppositely charged 5- 
bead surfactant molecules with one charged head bead 
and four neutral possibly hydrophobic tail beads. For 
finite values of the Coulombic and hydrophobic interac- 
tions we observe in general one or more aggregates or 
clusters of molecules. To analyze these structures in de- 
tail we define the clusters by assuming that any two beads 
belong to the same cluster if their distance is less than 
2(7. Among these clusters we identify the complex as the 
one that contains the polyelectrolyte chain. We moni- 
tor the number of clusters and the sizes of the complex 
and of the largest aggregates that are detached from the 
polyelectrolyte. These latter aggregates we identify with 
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FIG. 6: (a) The total mass (number of beads) of the com- 
plex and (b) of the largest detached micelle as function of the 
Coulomb strength Ab for different values of the hydropho- 
bicity parameter e L j = 0(D), 1(B), 2(o), 4(»), 6(A), 8(a) from 
above. 



isolated micelles. The number of identified clusters is 
shown in Fig. [5] as a function of the Coulomb interaction 
parameter Ab for different values of the hydrophobicity 
strength Elj- For vanishing interactions clustering oc- 
curs only by chance and almost each cluster observed 
represents one of the 33 single molecules in the system. 
For increasing interactions the molecules condense to a 
few larger clusters with only one or two remaining in the 
strong interaction case. When the Coulomb interaction 
vanishes, i.e. Ab = 0, the polyelectrolyte will always be 
detached from the micelles and thus a minimal number 
of two clusters remain even for strong hydrophobic at- 
traction. For any finite Coulomb interaction and high 
hydrophobicity on the other hand all molecules condense 
to a single complex in the present setup. The sizes of the 
complex and of the largest detached micelle are shown 
in Fig. where we plot the masses of these aggregates 
as function of the Coulomb parameter Ab for different 
values of the hydrophobicity Elj- It is clearly seen from 
these plots that a large detached micelle appears only for 
small Coulombic strength Ab < 1 and high hydrophobic 
attraction elj > 5. For these high elj > 4 the largest 
detached micelle shrinks rapidly with increasing Ab while 
at the same time the mass of the complex grows such that 



the surfactant is nearly quantitatively condensed in the 
complex for Ab > 1. For smaller hydrophobicity elj < 4 
on the other hand no larger detached micelles appear at 
all while the mass of the complex grows only moderately 
for increasing Coulombic strength Ab- Thus, while for 
high hydrophobicity the size of the largest detached mi- 
celle is inversely related to the size of the complex, the 
absence of large micelles for smaller hydrophobicity may 
be explained with their instability due to the Coulombic 
repulsion between the charged surfactants in combina- 
tion with entropy. As shown in the appendix the total 
binding energy of the surfactant tails in the micelle grows 
with proportionally to its mass M m j c with a surface cor- 

2/3 <• 

rection proportional to M m ' ic . The electrostatic energy of 
the surfactant heads on the other hand is proportional to 
^mlc wnne the entropy term behaves like M mic lnM m i c . 
The qualitative behavior of these three terms together 
with the approximate prefactors derived in the appendix 
may explain that for the finite number of available sur- 
factants the formation of large detached micelles is unfa- 
vorable for €lj < 6 and even for €lj = 6,8 the mass of 
the largest detached micelle is significantly smaller than 
the total mass M™!?f = 160 of available surfactant. 



C. Shapes of the complex 

In the case of complex formation we characterize the 
shape of the complex in terms of its matrix of inertia 



(a) (P)\ _ ~Yl 

r) i a ) = r 2 d a p 



(5) 



jec 



where the sum is over all M beads j that are part of 
the complex and Tj — (r^ , r^ 2 \ r^) is the position of 
bead j with respect to the center of mass of the com- 
plex. Diagonalizing Q( g results in the three eigenvalues 
fix > Ay > $z with respective axes x, y, z of the principal 
moments of inertia with the center of mass as the origin. 
Note that in these coordinates one has e.g. 



"d 7 = x 2 



(6) 



such that the radius Rg of gyration of the complex is 
given by R% = (& x +d y +'& z )/2 = ^TrO. The relative val- 
ues of the principal moments determine the global shape 
of the complex. The averaged ratios Si = (i?j//$x) and 
S2 = (^z/^x) are plotted as functions of the Coulomb 
strength Ab for different hydrophobicities Elj in Fig'- El 
These serve as indicators for the transition of the com- 
plex from a cylindrical shape to an almost spherical one. 
The cylindrical shape is characterized by two large ap- 
proximately equal inertia, i.e. Si « 1 and a third smaller 
moment such that S2 rs 0.2. ..0.3. For an elongated cylin- 
drical shape the z-axis corresponding to the smallest mo- 
ment $z is the symmetry axis. The spherical shapes in 
turn are characterized by three almost equal moments 
of inertia with S2 ~ Si « 1. We observe the following 
- see Fig. {7\ At X B = 0,e L j = only the PE chain 
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FIG. 7: (a) The ratio Si = fty/fix of the second largest 
and largest principal moments of inertia and (b) the ratio 
E2 = fix/fix °f the smallest and the largest principal mo- 
ments of inertia as function of the Coulomb strength Xb 
for different values of the hydrophobicity parameter slj — 
0(D), !(■), 2(o), 4(», 6(A), 8(A) from below. 



FIG. 8: The radius of gyration of the complex as function of 
the Coulomb strength A b for different values of the hydropho- 
bicity parameter e L j = 0(D), 1(H), 2(o), 4(«), 6(a), 8(a) from 
above. 



intermediate states in more detail below. A similar divi- 
sion into strong, weak and intermediate hydrophobicity 
we find for the behavior of the square radius of gyration 
Rq which is plotted as function of Xb in Fig. [HJ For small 
hydrophobicity El j < 3 the behavior of Rq parallels that 
of the square end-to-end distance R 2 as described above 
for the case Elj — 0. For high hydrophobicity elj > 5 
the radius of gyration of the complex shrinks already for 
small Xb ~ 1 to values corresponding to compact mi- 
cellar complexes. Again, the intermediate case Elj = 4 
behaves special: For small Coulomb strengths Xb < 3 the 
radius grows due to the stretching of the polyelectrolyte 
chain and only for higher A^ > 6 it again shrinks below 
the polymer coil value but still remains well separated 
from the high hydrophobicity case. 



is part of the 'complex'. The corresponding values we 
find for the neutral chain are Si = {fly/fix) = 0.87 and 
£2 = {$z/$x} — 0.34. These may be compared with 
{& y )/{& x ) = 0.87 and (0*)/{i? x ) = 0.27 calculated from 
data found for self avoiding walks j6^. Note that the 
relation {•& z / , d x ) > ($z)/($x) is in accordance with the 
Schwartz inequality for averages. For weak hydropho- 
bicity Elj < 3 the shape of the complex remains cylin- 
drical even for high Coulomb strength, while for strong 
hydrophobicity Elj > 5 the transition from cylindri- 
cal to spherical shape occurs already for small Coulomb 
strengths Xb ~ 1. For the intermediate hydrophobicity 
Xb = 4 though, we observe a competition between the 
formation of the cylindrical bottle brush and the spheri- 
cal micellar complex. This fact is expressed by the non- 
monotonic behavior of the shape parameter £2 as a func- 
tion of Xb ■ While for small Xb the cylindrical shape wins 
and £2 decreases with respect to the value of the neu- 
tral situation, a transition occurs between Xb = 2 and 
Xb = 4 where the shape parameter indicates a transition 
towards the spherical micellar complex, with £2 increas- 
ing to 0.6 and further towards 0.8 for As > 8. To verify 
the transition we analyze the internal structure of these 



D. Internal structure of the complex 

While the overall shape of the complex is characterized 
by the inertia matrix, its internal structure can be ex- 
plored by either the internal monomer- monomer correla- 
tions or by the density distributions of the different com- 
plexes in a given reference frame. We have used both ap- 
proaches to identify the key features of the self-assembled 
structures. To verify the idea of the molecular bottle 
brush like structures in the case of low hydrophobicity 
Elj and high Coulomb strength Xb, we have measured 
in particular the correlation between the polyelectrolyte 
monomers and the surfactant tail ends. As shown in Fig 
. using also data of ref. 0] we indeed find that the 
surfactant tail end is found in this situation with highest 
probability at a distance from the polyelectrolyte chain 
that increases with the length of the tail in accordance 
with the bottle brush picture. To monitor, on the other 
hand, the internal reorganization along with the tran- 
sition from the cylindrical bottle-brush to the spherical 
micellar aggregate, we chose as a reference frame a coor- 
dinate system defined by the directions of the principle 




FIG. 9: The pair correlation function g(r) for pairs of a 
polyelectrolyte monomer and a surfactant tail monomer at 
Bjerrum length \b = 8a. The functions are scaled to a max- 
imum value of 1. The maximal number of tail ends is found 
at r = 6cr from the polyelectrolyte backbone in the case of 
n = 5(U) monomer surfactants. For n — 10(a) this most 
probable distance is r = 9<r. 



moments of inertia of the complex with the origin at its 
center of mass. We define the x, y, z-directions such that 
the corresponding moments are ordered by t d x >'dy>'dz. 
As mentioned above, the z-axis then is the symmetry axis 
for the averaged cylindrical structures. 

In Figs. I1UI12I we show exemplarily for typical con- 
formations of the complex the density distributions mea- 
sured in this coordinate system. For vanishing hydropho- 
bicity £lj = and Xb — 4 the corresponding distribu- 
tions are shown in Figs. The densities as function of 
the distance r xy from the z-axis as plotted in Fig. HOf a) 
show clearly that the polyelectrolyte monomers and the 
surfactant heads are localized near the symmetry (z)-axis 
while the density of the surfactant-tail monomers sur- 
passes that of the former two species at larger distances 
from the central axis. The distributions measured along 
the z-axis on the other hand show the remarkable feature 
of maxima near the ends of the cylinder. Surprising as it 
may seem, this effect is even exhibited by the equilibrium 
coil of a neutral polymer chain and in particular in our 
simulation for vanishing interactions. 

The spherical micellar complexes, on the other hand, 
which we observe for sufficiently large parameter values 
of As and Elj are characterized by density distributions 
along and perpendicular to the z-axis that are nearly 
identical - see Fig. ^2 Only the different differential vol- 
umes used for measurement lead to deviations. A more 
natural variable in the spherical case is of course the dis- 
tance r from the center of mass. The density p(r) as 
displayed in Fig. El parallels the corresponding results 
of Wallin and Linse obtained by a self consistent spher- 
ical lattice field approach ^7j : The core of the micelle 
is formed by almost densely packed surfactant tails sur- 
rounded by surfactant heads that stick out of the surface 
and another layer that contains the polyelectrolyte chain. 

A transitional structure found for the values Xb = 4 
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FIG. 10: Density distributions of monomers in directions 
(a) along the z-axis and (b) perpendicular to the z-axis for 
Xb — 0, elj — 4 for polyelectrolyte beads (•), surfactant 
heads (o) and tail monomers (□). 



and Elj = 4 of the interactions near to the boundary be- 
tween the two states is investigated in Fig. ^| Although 
the structure is much more loosely bound it rather paral- 
lels the micellar aggregate displaying however a less pro- 
nounced layering in the spherical density distribution, see 
Fig. The comparison of Figs. El ( a ) an d (b) how- 
ever shows that the distributions display deviations from 
spherical symmetry. Nonetheless, the present data indi- 
cates that the transition between the two states of the 
complex is rather sharp. 



V. CONCLUSIONS 

In extensive MC simulations we have investigated the 
behavior of a single polyelectrolyte chain that interacts 
with oppositely charged surfactant molecules taking into 
account both long range Coulomb and short range hy- 
drophobic interactions. Exploring the parameter space 
in both of these interactions we have established a "phase 
diagram" for this system and several distinct regions 
( "phases" ) with characteristic structural behavior of the 
system as sketched in Fig. ^ f° r weak interactions both 
constituents behave nearly independently and the PE 
chain has a possibly stretched coil conformation. For high 
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FIG. 11: Density distributions of monomers in directions (a) 
along the z-axis, (b) perpendicular to the z-axis, and (c) as 
function of the distance from the center of mass for As = 16, 
el.j = 8 for polyelectrolyte beads (•), surfactant heads (o) 
and tail monomers (□). 
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FIG. 12: Density distributions of monomers in directions (a) 
along the z-axis, (b) perpendicular to the z-axis, and (c) as 
function of the distance from the center of mass for As = 4, 
el.j = 4 for polyelectrolyte beads (•), surfactant heads (o) 
and tail monomers (□). 



Coulomb strength but weak hydrophobicity we observe 
the formation of a complex with a cylindrical bottle brush 
structure in which the PE chain is stretched along the 
cylinder axis. If both interactions are strong, an inverted 
spherical micellar complex structure emerges with the SF 
tails constituting the core of a spherical micelle and the 
SF heads and the PE chain confined to the surface of this 
sphere. For intermediate values of the hydrophobicity 
there is a competition between the two structurally dif- 
ferent complexation modes resulting in a non-trivial be- 



havior of the shape as function of the Coulomb strength. 
For high hydrophobicity, but low Coulomb interaction 
the surfactant forms micelles to which the PE chain is 
only loosely attached. As we have reported in Ref.^fJ 
complexation in this system is in general weakened by the 
addition of salt. Our results can in principle be verified in 
experiments by systematically varying the surfactant tail 
hydrophobicity and the surfactant/polyelectrolyte charge 
and/or the dielectric permittivity of the solvent. 

A number of extensions of our present approach are 
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possible and would be interesting to follow. Of partic- 
ular interest would be a study of a multi PE-chain sys- 
tem. For small numbers of chains a corresponding simu- 
lation appears to be feasible using our present approach. 
One might expect to observe in such a case a competi- 
tion between the cooperative assembly of a multichain 
complex and single chain complexation. At high densi- 
ties in systems with many chains one expects the forma- 
tion of layered structures with many different possible 
phases 70] . However, high density and a large number 
of chains are out of reach of the present approach, due 
to the requirement of calculating the long range interac- 
tion between any pair of charged beads and the slowing 
down of the simulation by the caging effect at high den- 
sities. To simulate such systems either one has to release 
the hard sphere non-penetrability |7l| or resort to lattice 
methods and simulate only with a screened shorter range 
Coulomb interaction. 

Other situations that can lead to significant changes 
in the behavior might be found for a choice of differ- 
ent valency of the SF heads and the PE chains to study 
over-charging effects of the polyelectrolyte for high va- 
lency salt ions |72| . Also, in general, introducing neutral 
beads to the PE chain and varying the distribution of 
charges along the chain may lead to significant effects. 
In addition, the chain stiffness has been shown to be an 
important parameter. For a single micelle modeled as a 
charged sphere, many recent studies have revealed dif- 
ferent topologies of the defect de pen din g on the chain 

stiffnessjp M 0|3l m M El 00 

00mr 
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VII. APPENDIX 

A crude argument for the stability of micelles built 
from our model surfactants can be derived from the fol- 
lowing considerations: We assume that the core of a 
spheric micelle is given by mt tail monomers at a packing 
fraction r\ such that the core radius rc is given by 

.•.<>».> = f (^) " 3 . en 



The number of beads in the surface layer is then given 
by 

m a (m t ) = [r 3 c {m t ) - (r C (m t ) - ^f} V . (8) 

Assuming that each bead inside the core has 12 neighbors 
at an optimal distance (i.e. at the minimum of the poten- 
tial Vlj), while the surface beads have only 6 neighbors 
we find the total Lennard- Jones energy of the micelle as 

1 3 
——E LJ (m t ) = -(m t - m a (m t )/2)-£ LJ . (9) 
k B l 2 

The total Coulomb energy of the SF heads distributed 
on the micellar surface may be derived from the numer- 
ical solution of the Thomson problem for distributing k 
charges on a spherical surface which is given by the fit 
formula [73 

^=E c (k,m t ) = -^-^(l-ak-^+bk-^), (10) 
ksi rc(m t ) 2 

with a = 1.10461 and b = 0.137. Finally, we may es- 
timate the relevant entropy term by assuming that for 
every molecule it is given by restricting it to the volume 
of the micelle. This volume entropy term is thus 

^-S(k,m t ) = k\n(r(m t f /tt) (11) 
ks 

where Q is the total system volume. For our system of 
32 surfactant molecules with one charged bead and four 
tail beads we have mt = 4fc and a minimum of the free 
energy F(k,mt) with 

-^-F(fc, mt) = E c (k, mt) + E LJ (m t ) - S(k, m t ) (12) 
k B T 

at n ^ exists for a given Coulomb interaction only 
for sufficiently large hydrophobicity. Note, however that 
our estimate is too crude to give more than an order of 
magnitude result showing that Elj should be at least of 
the order of Xb and of the order of Elj ~ 1 for Xb = 
to find micellar aggregation for our case of surfactant 
length n — 5, i.e. m t = 4fc and the given overall number 
of surfactant molecules k < 32. 
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